"Dark" Z implications for Parity Violation, Rare Meson Decays, and Higgs Physics 
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General consequences of mass mixing between the ordinary Z boson and a relatively light Zd 
boson, the "dark" Z, arising from a U(l)d gauge symmetry, associated with a hidden sector such 
as dark matter, are examined. New effects beyond kinetic mixing are emphasized. Z-Zd mixing 
introduces a new source of low energy parity violation well explored by possible future atomic 
parity violation and planned polarized electron scattering experiments. Rare K(B) meson decays 
into n(K)£ + £~ (£ — e, /i) and -k{K)vv are found to already place tight constraints on the size of 
Z-Zd mixing. Those sensitivities can be further improved with future dedicated searches at K and 
B factories as well as binned studies of existing data. Z-Zd mixing can also lead to the Higgs decay 
H — > ZZd, followed by Z — > £^£\ and Zd — > i^^i or "missing energy", providing a potential hidden 
sector discovery channel at the LHC. An illustrative realization of these effects in a 2 Higgs doublet 
model is presented. 



I. INTRODUCTION 

The existence of cosmic dark matter is now essentially 
established. It appears to constitute about 22% of the 
energy-matter budget of the Universe, significantly more 
than the 4% attributed to visible matter Q|. Neverthe- 
less, the exact nature of dark matter remains mysteri- 
ous. Is it mainly a new, cosmologically stable, elementary 
particle that interacts with our visible world primarily 
through gravity or does it have weak interaction proper- 
ties that allow it to be detected at high energy accelera- 
tors or in sensitive underground cryogenic experiments? 
Both avenues of exploration are currently in progress. A 
discovery would revolutionize our view of the Universe 
and the field of elementary particle physics. 

Recently, a possible generic new property of dark mat- 
ter has been postulated [|2| to help explain various astro- 
physical observations of positron excesses M. The ba- 
sic idea is to introduce a new U(l)d gauge symmetry 
mediated by a relatively light Zd boson that couples to 
the "dark" charge of hidden sector states, an example of 
which is dark matter. Such a boson has been dubbed the 
"dark" photon, secluded or hidden boson, etc j| . Within 
the framework adopted in our work, however, we refer to 
it as the "dark" Z because of its close relationship to the 
ordinary Z of the Standard Model (SM) via Z-Zd mix- 
ing. Consequences of that mixing will be explored in this 
paper, where after describing the basic characteristics of 
the dark Z ', we provide constraints on its properties im- 
posed by low energy parity violating experiments such 
as atomic parity violation and polarized electron scat- 
tering. Future sensitivities are also discussed. We then 
briefly describe bounds on the mixing currently obtained 
from rare K and B decays along with the potential for 
future improvements. 

Perhaps the most novel prediction from Z-Zd mixing is 
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its implications for high energy experiments. In particu- 
lar, it leads to a potentially observable new type of Higgs 
decay, H — > ZZd, with pronounced discovery signatures 
that we describe J5[. We also discuss a 2 Higgs doublet 
(2HD) model that exhibits all the features of our general 
Z-Zd mixing scenario. (Some works of similar spirit, but 
different contexts can be found in, for example, Rcfs. |6| 
0-) 



II. SET UP 

We begin with what might be called the usual "dark" 
boson scenario. It is assumed that a new U{l)d gauge 
symmetry of the dark matter or any hidden sector in- 
teracts with the SU(3) C x SU(2) L x U(1) y of the SM 
via kinetic mixing between U(l)y and U(l)d EM. That 
effect is parametrized by a gauge invariant B liV Z d tv in- 
teraction 
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with e a dimcnsionless parameter that is unspecified (the 
normalization of the term proportional to e has been cho- 
sen to simplify the notation in the results that follow). 
At the level of our discussion, e is a potentially infinite 
counter term necessary for renormalization. Its finite 
renormalized value is to be determined by experiment. 
In most discussions, e is assumed to be < C(fcw x 10 -3 ). 
It could, of course, be much smaller |12[| . 

After removal of the e cross-term by field redefinitions 
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for the photon and Z boson fields, one is left with an 
induced coupling of the Zd to the usual electromagnetic 
current (with summation over all charged quarks and lcp- 
tons) 






(4) 



where the ellipsis includes W^ current terms and Qf 
is the electric charge (Q e = — 1). (It is generally as- 
sumed that U(l)d is broken and Zd becomes massive via 
a scalar Higgs singlet or a Stucckclbcrg mass generating 
mechanism |KJ, fhij.) Note also that the induced coupling 
of Zd to the weak neutral current via Eq. (||) is highly 
suppressed at low energies in the above basic scenario 
because of a cancellation between e dependent field re- 
definition and Z-Zd mass matrix diagonalization effects 
induced by e (see, for example, Ref. [^5J and our Appen- 
dices A and B). 

The phenomenology of the interaction in Eq. (m has 
been well examined as a function of mz d and e {e.g. 
Refs. @-@). With the assumption 10 MeV < m Zd < 
10 GeV and e < C(few x 10~ 3 ), bounds have been 
given and new experiments are underway to find the Zd 
via its production in high intensity electron scattering 
|jl9| . We will consider this same mass range for our phe- 
nomenological analysis in this work. The lower bound 
m z d ^ 10 MeV is required in order that astrophysical 
and beam-dump processes do not severely constrain the 
interactions of dark Z which, as discussed below, devel- 
ops an axionlike component for mz d — > 0. 

Because of its coupling to our particle world via the 
small electromagnetic current coupling in Eq. (0), Zd is 
often called the "dark" photon (even though that name 
was originally intended for a new weakly coupled long- 
range interaction p0|). 

Here, we generalize the above U{l)d kinetic mixing sce- 
nario to include Z-Zd mass mixing by introducing the 
2x2 mass matrix 
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where mz d and mz (with m 2 z <C m^) represent the 
"dark" Z and SM Z masses in the limit of no mixing. 
The Z-Zd mixing is parametrized by 
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with S a small model dependent quantity. We ignore the 
e contribution from Eq. (0) in the mass matrix, since its 
inclusion would affect this part of our discussion only at 
C(e 2 ) (see Appendix B). The assumed off-diagonal mz d 
dependence in Eq. (|6j) allows smooth mz d — > behavior 
for all e^-induced amplitudes involving Zd, even those 
stemming from nonconserved current interactions. Also, 
for simplicity, ordinary fermions are assumed to be neu- 
tral under U(l)d, ie. they do not carry any fundamen- 
tal dark charge. Their only couplings to Zd are induced 



through e and £z- More general cases are possible and 
interesting, but beyond the scope of this paper. 

So far, 8 is rather arbitrary, although < 5 2 < 1 is 
required to avoid an infinite-range or tachyonic Zd- One 
expects S to be small because of the disparity of mz and 
m z d ■ We later show that low energy phenomenology ac- 
tually requires S 2 < 0.006, while rare K and B decays 
have sensitivity to S 2 < 10~ 4 — 10~ 6 for low mass Zd- 
We will also demonstrate how the form in Eq. d|) natu- 
rally emerges in a simple 2HD extension of the SM, the 
details of which will be discussed in Appendix B. How- 
ever, we emphasize that our general results follow from 
Z-Zd mixing through a generic mass matrix of the form 
in Eq. (||) and are not exclusively tied to any specific ex- 
panded Higgs sector. That mixing could, for example, 
potentially arise from loop effects or dynamical symme- 
try breaking. 

Overall, mixing leads to mass cigenstatcs Z and Zd 
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where (see Appendix B) 
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It is expected that sin£ is very small (partly because of 
the assumed smallness of mz d /mz and partly because 
of small 6) and does not measurably affect Z pole pa- 
rameters (such as mz and T z) because these are shifted 
fractionally at 0(e z ), and require only £z < O(0. 01). 
However, it can, nevertheless, lead to other interesting 
new phenomenology which overcomes the mz d /mz sup- 
pression in £z- 

As the first example, we consider very low Q 2 parity vi- 
olating effects where the smallness of mz d /mz in the in- 
duced Zd couplings is offset by the m 2 z /m 2 z enhancement 
from Z vs Zd propagators. Then we describe the induced 
decays K — > irZd and B — > KZd, as well as the high en- 
ergy decay H — > ZZd, where the small induced coupling 
factor mz d /mz is overcome by ma /mz d , mB/mz d and 
m-H /mz d enhancements, respectively, in the longitudinal 
polarization component of the Zd production amplitudes. 



III. ATOMIC PARITY VIOLATION AND 
POLARIZED ELECTRON SCATTERING 

We begin our analysis by writing out the full Zd cou- 
pling to fermions from £ as well as £z- 
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where J^ m is given in Eq. (||) and 
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J" C = E( T 3/ - 2Qf sm 2 e w )f^f - T 3f h^ 5 f (10) 
/ 



with T 3f = ±1/2 (T 3e = -1/2) and sin 2 9 W ~ 0.23 is 
the weak mixing angle of the SM. The inclusion of Z- 
Zd mixing has introduced parity violation. The J^ Z% 
coupling is similar to the J^ Z^ coupling of the SM 
Z but reduced by ez in magnitude. Hence, the name 
"dark" Z, since it is the ez induced interactions that we 
primarily address. Note that the effects of e and ez can 
be combined into a simple form 



2cos0 w ^ Zd 
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by the replacement J™ (sin 9\y) 
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in Eq. (|l0|). In that format, one can judge the relative 
importance of s in low energy Zd phenomenology. It 
depends on the size of (e/e^)(cos 9w I sm9\y)- For e very 
small, it has little effect, but will be significant if e ~ Ez- 
The new source of parity violation in Eq. (0) or 
Eq. (O), is particularly important for experiments 
at Q < m Zd where the Zd propagator can pro- 
vide an enhancement owing to m 2 z <C rn z . The 
overall effect for parity violating amplitudes -M^c = 
(G F /2V2)F(sm 2 e w ) in the SM is (in leading order) to 
replace 
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sin 2 8\y — > Kd sin 2 9w 
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with |21 
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It is quite plausible that in a more complete theory, 
£ oc (mz d /mz)S = Ez- Then, the effects from kinetic 
mixing and Z-Zd mixing become similar in form and 
magnitude. Here, we allow e to remain a separate in- 
dependent parameter. 

Assuming no accidental cancellation between the pd 
and Kd in Eq. (H), Cesium atomic parity violation cur- 
rently provides the best low energy experimental con- 
straint on those parameters over the entire approximate 
range of interest (10 MeV < mz d < 10 GeV) since 
Q 2 -C m z . The nuclear weak charge measured in atomic 
parity violation (to lowest order in the SM) is given by 
Qw — — N + Z(l — 4 sin 2 0w) which when compared with 
experiment probes new physics. There is excellent agree- 
ment between the SM prediction for the weak charge of 



Cesium (including electroweak radiative corrections) |2 



(16) 



Q™a| 3 Cs) = -73.16(5) 
and the experimental value W\ Uj] 

Q w p (lfCs) = -73.16(35). (17) 

Based on the shift due to e, £z and 5 

Q™ ~> -73.16(1 + S 2 ) + 220— 6 2 cos6 w sindw, 

£z 

the above agreement then implies the following con- 
straints 



S 2 (i-1.27—) 
ez 



< 0.005 (la) 



(18) 



S 2 < 0.006 (one-sided 90% C.L.), for e < e z .(19) 

For £ ~ £z, the constraints on 8 2 become diluted and 
the possibility of cancellation occurs if one tunes £J £z — 
0.8. (We note that the fine tuning £J£z — 0.8 is similar 
to a relation employed in Ref. [[3J to try and reconcile 
what appears to be discrepancies in dark matter search 
scattering experiments on heavy nuclei. However, such a 
scenario is significantly constrained by the bounds on 6 
described below.) 

An independent constraint primarily applicable to Kd 
because of its relative insensitivity to pd comes from par- 
ity violating polarized electron-electron Moller scatter- 
ing asymmetries |28|, g9[. Experiment E158 at SLAC 
p(| measured the low energy value of sin 2 9 w(Q 2 ) at 
Q" ~ (0.16 GeV) 2 and compared it with expectations 
based on running the Z pole value sin dwijnz) down to 
low Q 2 (g9|. The good agreement with SM loop effects 
leads to (ignoring the small pd effect) 
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< 0.006 . 



(20) 



For m% ^$> (0.16 GeV) 2 and £z — £, the constraints 
in Eqs. (E9|) and (pFt) are essentially the same. How- 
ever, for a light m Zd < 200 MeV, the bound in Eq. @ 
can be somewhat diluted. Nevertheless, for some range 
of (s, mz d ) values, Eq. ( p0| ) can provide more restric- 
tive bounds on 5. For example, consider e ~ 2 x 10~ 3 
and mz d — 100 MeV which lie in the region favored by 
the current discrepancy between theory and experimen- 
tal values of the muon anomalous magnetic moment [pTJ . 
In that case, Eq. (EOT) becomes 



\S\ <0.01 



(21) 



which is considerably tighter than Eq. (|l9|). If the muon 
anomaly discrepancy is because of a light Zd and e ~ 
10~ 3 , that boson's effect on the value of sin 2 9w extracted 
from future more precise very low Q 2 parity violating 
experiments |32[ could eventually become observable. 





FIG. 1: Examples of diagrams contributing to b 
ilar diagrams give rise to s — > dZ d - 



sZ d . Sim- 



The sensitivity in Eqs. (^0|) and (^Tj) is expected to 
improve by up to an order of magnitude from ongoing 
and proposed polarized ep and ee scattering experiments 
at JLAB H as well as proposed Q 2 ~ (0.05 GeV) 2 ep 
studies at MESA in Mainz pjj. Our analysis illustrates 
the complementarity of direct searches at intense elec- 
tron scattering facilities in JLAB and Mainz for a light 
vector particle (the "dark" photon coupled through ki- 
netic mixing) produced via electron scattering, with low 
Q 2 measurements of sin 8\y hi parity violating experi- 
ments (that probe e and the mass mixing of the "dark" 
Z) . We also note that proposed measurements of atomic 
parity violation for ratios of different nuclear isotopes 
would eliminate atomic physics uncertainties as well as 
any dependence on pd |?4]-p7[ . They would then be sensi- 
tive to (e/ez)S 2 but with negligible Q 2 dependence (since 
Q 2 ~ 0). It is amusing to note that in principle, very low 
energy measurements of sin 6\y in atomic parity viola- 
tion and low Q 2 polarized electron scattering experiments 
could find different sin 9w results from one another if a 
very low mass Zd is contributing to both, because of the 
Q 2 dependence in Eq. (p"4|). 

Our conclusion, based on the above discussion, is that 
currently, 5 2 < 0.006 is a modest, reasonably reliable 
constraint for most values of rnz d , although hne tuning 
of e and ez could loosen the bound. That constraint can 
be much stronger for e ~ 10~ 3 [sec Eq. (J2l])], and could 
be further improved significantly by future low energy 
parity violating experiments. For now, the bound S 2 < 
0.006 provides a starting point for comparison with the 
sensitivity to 5 2 in rare K and B decays which we next 
describe. 



IV. RARE K AND B DECAYS 

Experimental studies of rare flavor-changing weak neu- 
tral current decays of K and B mesons have proven to be 
powerful probes of high and low scale "new" physics phe- 
nomena. Here, we illustrate the effect of Z-Zd mass mix- 
ing on the transition amplitudes s — > dZd and b — > sZd 
induced within the framework of Cabibbo-Kobayashi- 
Maskawa (CKM) charged current mixing (See Fig. [I]). 
Those loop induced couplings can lead to decays such as 
K —> irZd and B — > KZd or K*Zd characterized by the 
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FIG. 2: Z d lifetime with Z d mass for <5 2 = 10" 4 with e = 
(solid blue curve) and e = 2 x 10~ 3 (dashed blue curve) cases. 
We take p, <f>, J/ip, T masses as the representative threshold 
for decays to mesons. 



signature Zd — » £ + £~ (£ — e or p) with invariant mass 
mu = m z d or Zd — > missing energy where Zd decays 
into vv or essentially undetectable light hidden sector 
particles. In all such 2-body decays, the mono energetic 
outgoing 7T or K will provide a tight constraint (for a 
given mz d ) and a very distinct overall signal. 

Here, we note that the phenomenology of Zd is affected 
by its lifetime Tz d ■ A sufficiently large value of Tz d will 
allow Z d to escape the detector and lead to a missing 
energy signal. However, for smaller values of Tz d , a dis- 
placed vertex can provide a distinct signature. In Fig. 0, 
using representative values of 5 and e, we have plotted 
rz d for 10 MeV < mz d < 10 GeV, assuming that Zd 
only decays into SM final states. We provide a simple 
formula for the partial width of Zd into SM fermions, 
T(Z d -> //), in Appendix C. 

Of course, the amplitudes for dsZd and sbZd being 
loop induced will in general depend on the details of 
the complete model considered, including its underlying 
Higgs flavor symmetry breaking structure. Those details 
are beyond the scope of this paper where we are pri- 
marily interested in the generic effects of Z-Zd mixing 
parametrized by ez = { r mz d /mz)8 in Eq. (g). 

A simple illustrative example of a scenario that leads 
to Z-Zd mixing and CKM induced flavor-changing weak 
neutral currents is the Type-I 2HD model discussed in 
Sec. VI and detailed in Appendix B. There, the under- 
lying U(l)d gauge symmetry naturally forbids tree level 
flavor-changing neutral currents in the scalar and pseu- 
doscalar Higgs sectors. It also yields, through Higgs dou- 
blet and singlet vacuum expectation values, a mechanism 
to provide mass for Zd and give rise to a small <5 in Eq. (|g). 

To obtain the induced Zd flavor-changing amplitudes, 
we can make use of existing CKM loop induced cal- 
culations for dLjfj.SLZ 11 and s^^bhZ^ amplitudes |3S|1 
and replace Z —> ezZd- (See Fig. |jj) (We ignore ki- 
netic mixing induced couplings, since their effects are 
highly suppressed. For example, Rcf. |}9l found BK(B — > 



KZd) ~6x 10~ 7 e 2 for mz d — 1 GcV. As we demon- 
strate, mass mixing, sz, induced rates can be much larger 
and potentially observable.) As an alternative compu- 
tational strategy, if we are primarily interested in rel- 
atively light Zd bosons compared to rriK and tub, we 
can employ the Goldstone boson equivalence theorem 
j40| to obtain amplitudes for longitudinally polarized Zd 
bosons from flavor-changing axionlikc pscudoscalar cou- 
plings well documented in the literature. For our pur- 
pose, the latter approach will suffice; however, the direct 
Z calculation provides a nice cross-check. Nevertheless, 
we note that the results discussed below should be viewed 
as somewhat incomplete and should be taken as approx- 
imate. 

The relevant di7 M Si3' i a and SLlpbhd^a axion cou- 
plings were computed for the 2HD model more than 
30 years ago by Hall and Wise M and independently 
by Frere, Vermaseren and Gavela|42|. More recently, 
they were checked and applied to the decay B —} Ka. 
a — > £ + £~ in Rcf. J43]. Here, we use those results to 
estimate the branching ratios for K — > TiZd (longitu- 
dinal) and B — > KZd (longitudinal) which should ap- 
proximate the full Zd final state rates up to corrections 
of 0{m Zd /m 2 K ) and 0(m z /m/g) respectively. Compar- 
ison of those estimates with experiments can then be 
used to constrain S for the ranges m z -C (rriK ~ fn^) 2 
and Tnr z <C (ms — mx) 2 modulo regions not covered 
because of experimental acceptance cuts on the data 
(which arc beyond the scope of this paper). For example, 
mz d < 140 McV is not covered because of 7r° — > e + e~j 
Dalitz decay background. Similarly, masses of Zd near 
charmonium resonance regions are not covered. 

We begin with the predicted branching ratio for K — > 
irZd (longitudinal) in the 2HD model. Based on the anal- 
ysis in Ref. GIL but adjusting for a modern mt value, 
since top now dominates the amplitudes in Fig. n\ 
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(22) 



where the numerical factor in that expression includes 
QCD suppression effects and depends on the physical 
charged scalar Higgs mass of the 2HD model. Those un- 
certainties should be considered part of the overall model 
dependence of our analysis. 

The Zd produced in Eq. (22h is expected to decay 
promptly (see, however, Fig. 2|) to £ + £~ pairs with in- 
variant mass mz d or to missing energy that might be vv 
or light hidden sector particles. Those decays would add 
to the SM predictions and should be part of the experi- 
mentally measured branching ratios pL H3, Ej| 



BR(K + 
BR{K + 
BR(K + 



n + e + e 



cxp = (3.00 ± 0.09) x 10~ 7 (23) 
7rV + AOc X p = (9.4 ±0.6) x 10~ 8 (24) 
7T + ^)cx P = (1.7 ±1.1) x 10~ 10 (25) 



unless eliminated by acceptance cuts which would negate 
bounds in certain mz d regions. For example, the result in 
Eq. (p3) applied a m ee > 140 MeV cut while Eq. ([25]) was 



obtained with a rather stringent cut on E v . Clearly, a 
new round of bump hunting in the £ + £~ spectrum is war- 
ranted. Toward that end, we note that Zd — > l + (r de- 
cays will have a characteristic polarized spin-1 sin 9 dis- 
tribution relative to the longitudinal polarization of the 
Zd- Unlike the spin-0 axion case, where because of chiral 
conservation the a preferentially decays to the heaviest 
fermion possible and the distribution is isotropic, we ex- 
pect BK(Zd —> e + e~) ~ BR(Zd — » /.i + /j,~) modulo phase 
space. 

With the above caveats, we compare Eq. ( p2| ) with 
(P3|), (p4j) , and (|25| ) which agree with SM expectations 
and find rather tight bounds 



l<*| < 0.01/VBR(Z d -> e+e-) (26) 

\S\ < 0.001/v/BR(Z d -> missing energy) (27) 

modulo acceptance cut criteria. 
Eqs. © and (|J) yield g§ 



BR(Z d 
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BR(Z d -*■ vv) 



(28) 



where e from kinetic mixing now comes into play. For 
e ~^> Sz , the charged lepton decays dominate and Eq. (l2£ 



is more applicable. For e < Ez, the tighter constraint 
in Eq. ( p7| ) takes precedence. Of course, both should 
be used cautiously, given their model and experimental 
acceptance dependence. 

For the case of B — > KZd (longitudinal), we can apply 
a similar approach and find pl|-j43|| 



BR(B -»■ KZ d )ion g s 0-IS 2 



(29) 



The relatively large coefficient in Eq. (E9) results from a 
factor of mf in the b — > sZd loop induced correction from 
Fig. ||. That factor makes rare B decays a particularly 
sensitive probe of the Zd. Employing the recent bounds 
that follow from the discussion of B — > Ka, with the 
axion-type particle a — ¥ £ + £~ in Refs. |3!], Q implies 
conservatively BR(B -> KZ d -> K£+£~) < 10" 7 , while 
the bound from B-dccay containing missing energy are 
based on |, |t], |§ 



BR{B + -> K 
We then roughly find 



/cxp 



< 1.4 x 10 



-5 



(30) 



\S\< 0.001/ ^/BR(Z d 
\S\ <0M/^BR(Z d - 



+ £+£-) (31) 

missing energy). (32) 



It has been suggested |35| that even tighter bounds 
may be obtained from dedicated searches for £ + l~ pairs 
in B decays, particularly if displaced vertices result from 
suppressed decay rates. Nevertheless, even the relatively 
crude bounds in Eqs. ( pi] ) and (p3) are very constrain- 
ing where applicable and are likely to be significantly 
improved by future dedicated searches. 



On the basis of our analysis, it is clear that rare K and 
B decays provide sensitive windows to Z-Zd mass mixing 
and should be further explored in future high intensity 
experiments. In fact for both cases, a more refined binned 
analysis of existing data would likely result in tighter 
bounds than those in Eqs. (|26| ) and (|3l]) or even uncover 
a hint of the ZdS presence. Although applicable to a 
limited range of mz d and dependent on the Z4 branch- 
ing ratios, one can easily conclude \5\ < 0.01 — 0.001 over 
some restricted mz d domain. In addition, further im- 
provements are possible and warranted. That constraint 
on 8 sets a standard for other rare decay studies. As we 
show in the next section, it is possible that searches for 
the rare Higgs decay H — > ZZd have the statistical sig- 
nificance to also explore \8\ < 0.01 — 0.001 but have the 
potential advantage of covering a much broader range of 
mz d values including mz d > 5 GeV if backgrounds can 
be controlled. 



V. HIGGS DECAYS 

We now address a primary consequence of our paper, 
the decay H — > ZZd induced by Z-Zd mass matrix mix- 
ing. To put our analysis into a current day perspective, 
we take m# = 125 GeV, a value roughly suggested by 
early small excesses at the Large Hadron Collider (LHC) 
in the expected decay modes H — s- 77, WW*, ZZ* 
P9| , |50|| . We note, however, that our findings regard- 
ing the sensitivity of Higgs searches for H — > ZZd are 
fairly independent of the exact value of tuh ■ 

To set the stage, we estimate that, roughly, one ex- 
pects each LHC experiment to have about 75000 Higgs 
bosons in the existing data before cuts (for the integrated 
luminosity of 4.7 - 4.9 fb -1 with E c . m . = 7 TeV) for 
m H = 125 GeV in the SM. In Table |, we iist the expected 
Higgs decay branching ratios within the context of the 
SM. Of particular interest for comparison with H — > ZZd 
arc the SM decays (1) H -)■ ZZ* ->■ l$£[l$q and (2) 
H -> ZZ* -> £ + £~vv where the * signifies a "virtual," 
off mass shell boson and I = e, fi. The first of these, even 



H Decay Channel 


Branching Ratio 


bb 


0.578 


WW* 


0.215 


99 


0.086 


r+r- 


0.063 


cc 


0.029 


ZZ* 


0.026 


77 


2.3 x 10" 3 


Z 1 


1.5 x 10" 3 


H ->■ ZZ* -*■ £+£~£+q 


1.2 x 10~ 4 


H ->■ ZZ* -+l + £~vv 


3.6 x 10~ 4 



at the BR ~ 10~ 4 level, may have already been seen at 
the LHC where a handful of candidate events have been 
reported. If it truly is a Higgs signal, hundreds more 4- 
lcpton if £^£^£2 events will be clearly observed in the 
coming years. The second decay, H — > £ + £~vD, is more 
difficult and to our knowledge has not been experimen- 
tally studied. 

For the first case, one lepton pair will have an invariant 
mass of mz — 91 GeV while the second pair will have an 
invariant mass ranging from to about 34 GeV with a 
differential decay rate distribution as depicted in Fig. g. 
The second mode H -> ZZ* -> £ + £~vv, with the neu- 
trinos identified by missing energy, while experimentally 
more challenging should be searched for as well, since it 
can be used to constrain potentially invisible decays of 
the Zd, as we subsequently discuss. 

As we shall see, the decays H — > ZZd are significantly 
enhanced beyond naive expectations, even for very small 
mixing. To appreciate that phenomenon, we remind the 



reader that for a very heavy Higgs (m H ^> 



*W> 



r) the 



decay rates for H — > W + W~ and H — > ZZ can become 
enormous, growing like ~ g 2 m? H /my, V = W, Z with 
increasing m#. That behavior comes about because the 
final state W and Z bosons are longitudinally polarized, 
resulting in a ~ m 2 H /rriy enhancement factor at the de- 
cay rate level (for each final state gauge boson) . 

Such an effect is a manifestation of the Goldstone bo- 
son equivalence theorem which states that at high ener- 
gies (s ^> my), S'-matrix elements involving W ± and Z 
bosons arc equivalent, up to 0(my /\fs), to the corre- 
sponding amplitudes in the Higgs-Goldstonc scalar the- 
ory with the Goldstone boson replacing W/~ , Zl (longitu- 
dinal components). In the heavy Higgs limit, the W + W~ 
and ZZ decay products are essentially longitudinally po- 
larized and behave like their Goldstone boson compo- 
nents. The Higgs coupling to Goldstone bosons is of the 
form {—ig/2)m 2 H /mv, and squaring that coupling and 
dividing by 1/ttih gives the T(H -» VV) ~ g 2 m z H /m v 
exhibited by heavy Higgs decays. We note that the lon- 
gitudinal polarization of the gauge bosons can be very 
helpful in identifying a Higgs decay since the subsequent 
decay W or Z — >• leptons have a characteristic angular 
distribution oc sin 9 relative to the polarization. 

Of course, our example of 125 GeV Higgs is too light 
to decay into W + W~ or ZZ pairs. It can, however, 
decay into one real and one virtual boson with the latter 
directly producing a lepton pair with an invariant mass 
distribution as illustrated in Fig. |52| . The integrated 
partial width for H -» ZZ* -t tjiyi^l^ is, however, 
suppressed by a/47r (from the Z*£\£^ coupling and 3- 
body phase space) and the small BR(Z — > £ + £~) ~ 2 x 
0.034 for £ = e, /.1. One finds 



T(H -> ZZ* -+ £\£'{£\€ 1 ) ~ 1.8 x 10" 



G, 



8V2tt 



TABLE I: Standard Model Higgs decay branching ratios for 
m H = 125 GeV (T„ ~ 4.1 MeV) from Ref. §1. 



m H m w 

(33) 

with no significant sign of enhancement for longitudinal 
polarization, which is not surprising, since nitf/raz — 1.4 
in our example. Nevertheless, even with the 10 -6 sup- 



pression factor in Eq. (|33|), it is expected that a SM 
125 GcV Higgs should be starting to be seen with about 
several events per experiment in existing data, after ac- 
ceptance cuts, and with hundreds more to follow in sub- 
sequent years. So, Eq. ([33]) represents a decay rate stan- 
dard that is easily discernible if backgrounds are in check. 
We note that the decay rate for H -> ZZ* -> i + l~vi> 
is expected in the SM to be about 3 times larger than 
Eq. ( p3| ) but more difficult to measure. 

Now we come to the decay H — > ZZd owing to Z- 
Z c \ mixing in our "dark" Z scenario. That mixing, 
parametrized by ez = (mza/mz)^ a very small quan- 
tity, might naively appear to be negligible since it leads 
to a tiny HZZd coupling ~ (g/ cos6w)ir>-z£z- Conse- 
quently, the H — > ZZd decay rate will be suppressed by 
e 2 z = (mz d /mz) 2 5 2 . However, because of the Goldstonc 
boson equivalence theorem, we gain an enhancement fac- 
tor of ~ (rriH /iTiz d ) 2 in the decay rate for longitudinally 
polarized Zd final states (a feature that may also help in 
identifying their subsequent Zd — > i^2 products via an- 
gular distribution if statistics suffice). That enhancement 
negates the small mz d /vnz factor in the HZZd coupling. 
Also, there is no a/ An suppression for H — > ZZd-, only 
the small BR{Z ->■ £+£') ~ 2 x 0.034 that needs to be 
included for Z identification. A detailed calculation (see 
Appendix B) leads to 



T(H -> ZZ d -> till It t%) 
±7x 10~ 3 - 



1^* BR(Zd 



"2 *2 



(34) 



Note the rv? H behavior that results from Z and Zd be- 
ing produced in their longitudinal polarization modes. 
A similar formula with BR(Zd — > ^2^-2) replaced by 
BR(Zd — > missing energy) applies to the case Zd — > vv 
or invisible "dark" particles. 

In terms of its branching fraction relative to the SM 
expected width, one finds 



T(H -> ZZ d ) 



16 x S 2 < 0.1 



(35) 



r| M (125 GeV) 

with T| M (125 GcV) ~ 4.1 x 10~ 3 GcV §| and using the 
low energy bound in Eq. (|j~9|). We see that as much as 
10% of all LHC Higgs decays could be producing ZZd- 
With current statistics, even a 10% loss of SM expecta- 
tions would not be noticed; but eventually it would be 
uncovered by precision Higgs production and decay stud- 
ies. 



(36) 



Taking the ratio of Eqs. ( 
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and 
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FIG. 3: Differential decay rate H ->■ ZZ* -> Z£ + T -> U vs 
£ + l~ invariant mass with run = 125 GeV in the SM (in blue). 
For the illustration, H -* ZZ d -> Z£ + l~ with mz d = 5 GeV 
and <5 2 BR(Z d -> £ + £~) = 1(T 5 (which would need V H i ggs ~ 
10 6 for 3<r evidence) is also shown (spike at the 5 GeV bin in 
red). Bin size is selected to be 2 GeV. 



for invisible Zd decays. Even for 5 2 ~ 10 4 , well below 
the atomic parity violation bound of 0.006 in Eq. (|lS|), 



-1 <-l *-2 *-i 



or 



one would expect H — > ZZd events with 
£ + £~ + missing energy to be starting to appear or already 
present in LHC data. If there are no Zd — > dark parti- 
cles decays, we expect the branching fractions of Zd into 
£ + £~ to be given by Eq. (28). Therefore, in that case, one 
expects BR(Zd —> £ + £~) to be relatively large, particu- 
larly if (e/ez) 2 ^ 1. If Zd — >• dark particles dominates 
its decay rate and significantly dilutes BR(Zd — > £ + £~), 
one still has the possibility of seeing H — > ZZd — > £ + £~ + 
missing energy, although this perhaps is more experimen- 
tally challenging. Of course, given the original motiva- 
tion for introducing a Zd into astrophysics as a way of 
explaining positron excesses through its decays, a rela- 
tively large BR(Zd — > £ + £~) might be expected. 

Returning to Eq. (|3g), we see that even for a some- 
what suppressed BR(Zd — > £ + £r), the LHC experiments 
should be able to search for a Zd in the H — > t^t^H^t^ 
decay chain down to S 2 BR(Z d ->• £ + £~) ~ O(10" 5 ), de- 
pending on backgrounds. (The domain explored by rare 
K and B decays for some subset of mz d values.) The 
signature, two isolated lepton pairs £f£^ + l^t^ with 
a total invariant mass of tuh and individual masses of 
mz and mz d should stick out as a spike in the invari- 
ant mass plot of Fig. pi as illustrated for mz d = 5 GeV 
and 5 2 BR{Z d -> £ + £~) = 10~ 5 . In the bin centered 
at Mu — 5 GeV, the SM expectation from Higgs of 
m H = 125 GeV is - 6.3 x 10" 9 GeV, while the sig- 
nal associated with H -^ ZZd is ~ 4.5 x 10~ 8 GeV. 
With existing data of Njj — 75000, no meaningful num- 
ber of signal or background events are expected, and one 
would need Nmggs — 10 6 for 3a evidence (beyond the 
SM H ^ ZZ* -T A£ channel) at the LHC experiments. 
However, this simple estimate ignores other reducible and 
irreducible backgrounds and a more reliable statement re- 
quires inclusion of such details. Also, the £^£2 decay pair 



from Zd should exhibit an angular distribution consis- 
tent with its longitudinal polarization. That sensitivity 
is potentially orders of magnitude below the 8 2 < 0.006 
already established by atomic parity violation. We note 
that while the Higgs decay constraints on 8 may not sur- 
pass those derived before from rare K and B decays, they 
are applicable well beyond the O(GeV) regime of mz d , 
relevant for the meson decays. They represent a poten- 
tially unique broad capability of the LHC unmatched by 
low energy experiments. 

We should point out that current searches for H — >■ 
ZZ* —$■ 4£ are likely to miss H — > ZZd because they 
generally cut out a lighter second lepton pair with M« < 
15 GeV, i.e. the range of interest, in order to avoid Zj* 
backgrounds. Hopefully, our results will provide some 
incentive for revisiting the low mass region in search of 

z d . 

In addition to Zd — >• £ + £~ , one should mount a search 
for H — > ZZd -^ £ + £~ + missing energy. Here, one might 
be helped by the fact that the missing energy and miss- 
ing momentum of the Zd decay pair are nearly equal. A 
thorough study of LHC capabilities for uncovering that 
decay mode is clearly warranted. We also add that the 
Higgs can have a decay mode H — > ZdZd, in our frame- 
work. The rate for this decay is proportional to 8 4 , so, 
roughly, it is suppressed compared to the ZZd mode by 
0(S 2 ) which, given our bound in Eq. (|l9|), is a suppres- 
sion of 0.006 or smaller. The rate for the ZdZd channel 
could be enhanced if hidden sector scalars that couple 
directly to Zd and give it mass are allowed to mix with 
the SM sector Higgs scalars. 



VI. A 2 HIGGS DOUBLET EXAMPLE 

In the preceding discussion, we examined the dark Z 
phenomenology in a general framework. As mentioned 
before, the main ingredient we introduced was mass mix- 
ing between the SM Z and the Zd which could be real- 
ized in a variety of models. In this section, to demon- 
strate how our general framework might be realized, we 
will consider a 2 Higgs doublet extension of the SM. (See 
Ref. |53| for a recent review on 2HD models.) Here, we 
assume two SU(2)l x U(1)y Higgs doublets, Hi and H 2 , 
but allow H 2 to carry a "dark" charge that couples it di- 
rectly to U(l)d- Note that the assumption of the U(l)d 
in our example is well-motivated, as it allows the model 
to evade severe constraints from flavor-changing neutral 
currents that arc often addressed through the introduc- 
tion of a Z2 symmetry in generic 2HD models. We also 
allow, for generality, a singlet scalar, Hd, that also pro- 
vides part of the Zd mass through its "dark" sector vac- 
uum expectation value Vd- 

With the above assumptions, H 2 does not couple di- 
rectly to ordinary fermions, but does contribute to W^, 
Z and Zd masses as well as Z-Zd mixing through its vac- 
uum expectation value v 2 . Such a setup is akin to what 
is often called a Type-I 2HD model p4|. Here, we will 



take H\ to be a SM-likc Higgs scalar, identified as H 
in our preceding general analysis. To keep the discus- 
sion simple, we ignore scalar mixing among the Hi, H 2 , 
and Hd states. The V\, v 2 , and Vd vacuum expectation 
values of Hi (the SM doublet), H 2 and Hd give rise to 
8 = sin /3 sin /3d where tan/3 = v 2 /v\ and tan fid = v 2 /vd, 
as will be shown in Appendix B. The condition of a 
SM-like Hi can be satisfied, to a good approximation, 
for tan/3 < 1/3, and does not require a large hierar- 
chy of scales in the Higgs sector. The constraints on 8 
previously discussed will however constrain the product 
sin /3 sin (3d ■ 

There are many additional features of our 2HD model 
worth studying. For example, nonzero Higgs scalar mix- 
ing (which we set to zero) could give rise to enhancements 
in H — > ZdZd, as mentioned before, or perhaps H — > hh 
(h being a lighter Higgs scalar remnant of H2) | 55 . Those 



possibilities are interesting but more model dependent. 



VII. SUMMARY AND CONCLUSION 

In this work, we explored the possibility of mass mix- 
ing between the Z boson of the SM and a new light vector 
boson Zd associated with a hidden or dark sector U(\)d 
gauge symmetry. Such a light state has been invoked 
in discussions of astrophysical anomalies that may origi- 
nate from cosmic dark matter. We dub this new vector 
boson the "dark" Z, as its properties are analogous to 
that of the SM Z. In particular, the couplings of Zd can 
provide new sources of parity violation and measurably 
affect the decay of the Higgs through novel channels such 
as H — > ZZd- Existing atomic parity violation, polarized 
e scattering, and rare K and B decay data already place 
interesting bounds on the degree of Z-Zd mass mixing, 
but further improvement is possible and warranted (see 
Table H 1 . 

The presence of kinetic mixing affects the phenomenol- 
ogy of Zd, but much of the main physics discussed in our 
work persists even in the absence of kinetic mixing. Vari- 
ous experimental efforts arc currently devoted to possible 
signals of the "dark" photon, based solely on the possibil- 
ity of kinetic mixing between U(l)d and the SM photon. 
Here, we want to emphasize the m 2 z /m 2 z enhancement 
factor in low energy parity violation and the longitudi- 
nal polarization enhancement Ez d /mz d , with Ez d the 



1 One could contemplate searching for Z^ effects in precision 
neutrino neutral current cross section measurements such as 
v^e — > v^e or deep-inelastic u^N — > u^X. However, to be com- 
petitive with anticipated low energy parity violation polarized 
electron scattering or atomic experiments, those neutrino stud- 
ies would have to reach ~ ±0.1% statistical and normalization 
uncertainties, a challenging task that would likely require a high 
energy neutrino factory (see Ref. ]56|). A detailed discussion of 
Zd effects on neutrino cross sections will be given in a separate 
publication. 



Process 


Current (future) bound on S 


Comment 


Low Energy Parity Violation 

Rare K Decays 

Rare B Decays 

H^ZZ d 


\5\ < 0.08 - 0.01 (0.001) 

\S\ < 0.01 - 0.001 (0.0003) 

\S\ < 0.02 - 0.001 (0.0003) 

\5\ < (0.003 - 0.001) 


Fairly independent of mz d - Depends on e. 

m\ < m 2 z <S rn 2 K . Depends on BR(Z<j). 

m^ < m Zd <C m B . Depends on HK(Zd)- Some mass gap ~ 3 GeV. 

m % d ^ (ffiH — mz) 2 ■ Depends on BR(Z^) and background. 



TABLE II: Rough ranges of current (future) constraints on S from various processes examined along with commentary on 
applicability of the bounds. These processes have negligible sensitivity to pure kinetic mixing effects. 



energy of Zd, in rare meson decays and the Higgs de- 
cay H — > ZZd- These enhancements make such pro- 
cesses particularly sensitive to very small Z-Zd mixing. 
In particular, future polarized ep and ee scattering exper- 
iments can provide further probes of the scenario we have 
considered in this work. These parity violating probes 
are sensitive to a wide range of Zd masses, including 



m z d 



< 140 McV. where other searches fail because of 



7r Dalitz decays background and are independent of Zd 
branching fractions. The rare K and B decays currently 
provide some of the most stringent bounds on the degree 
of Z-Zd mixing, however they depend on the Zd branch- 
ing fractions and also do not apply to mz d above the me- 
son mass. In addition, there can be gaps in the bounds, 
for example in the mz d charmonium mass region. 

In the event of the discovery of a SM-like Higgs at 
the LHC, say at ~ 125 GeV based on current hints, a 
new front in the search for a dark Z can be established. 
The Higgs decay data are particularly unique for mz d > 



5 GeV, and hence probe a part of parameter space that 
is inaccessible to meson data. The reach for this new 
physics can be extended well beyond the current limits 
through precise measurements of Higgs decays, as may 
be done at an e + e~ or /i + /i~ collider if high statistics 
are available. We conclude that pushing the above types 
of experiments as far as possible is strongly motivated, 
for they could be windows to the "dark side" of particle 
physics. 
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Appendix A: Gauge Kinetic Terms 
The gauge kinetic terms allowed by the gauge symmetries SU(2)l x U(l)y x U(l)d are 
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1 



^aiiiy 



= — -B B^ v 

4 M 2 cos 0w 






1 



-Z. 
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rOftv 



dp.is d 



(Al) 



with F^ v = d^Fy — dyF^. The hatted quantities are fields before the diagonalization of the gauge kinetic terms. The 
diagonalization is done by the field redefinition known as a GL(2,R) rotation 



7° 
B„ 



^/l-e 2 /cos 2 9 w 
— e/cos6w 1 



7° 
B„ 



after which, B gets a Zd component proportional to e while Zd does not get any B component. 



— _ L R TZ^v _ _ 7° yOfiiy 
gauge — .-D/iiz-D A^djiv^d 



(A2) 



(A3) 



We will take Z® = Z® and B^ = B^ L + (ej cos 9w)Z^ and ignore 0(e 2 ) terms from here on. After electroweak 
mixing with Weinberg angle 9w 




we get 



cos t>w sin t>w 

sin 9w cos 9w 



A^ — A/j, sZ dft 

z f _ l = Z^+e tan 6wZ dil 




(A4) 



(A5) 



yd _ yO 



as an effect of the gauge kinetic mixing. Thus, Z° is unaffected to O(e) while both A^ and Z° are shifted by the 
gauge kinetic mixing followed by the electroweak mixing. However, the bare fields do not take into consideration 
Z°-Z d ) mixing via the mass matrix from the Higgs mechanism which we will deal with in the following. 



Appendix B: Scalar Kinetic Terms 



The scalar kinetic term is given by 



^scalar = /_^ \^^i\ 



(Bl) 



where i runs for all Higgs scalars. Considering only neutral components of gauge bosons, we have 

D^i = (fl M + ig'Y^B^ + igTz[§i}W^ + igaQafii]^) $< (B2) 

before gauge kinetic diagonalization where Y , T3, and Qd are hyperchargc, isospin, and dark charge, respectively. 
After symmetry breaking, the scalars can be written with the vacuum expectation values (vi). 

$ i = -^=(H i +v i ) (B3) 

1. Vector boson mass 



From Eq. (Bl), we can get the relevant vector boson mass terms 
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The mixing of two vector bosons is given by 



with 
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-m 2 7°7° 



cos£ — sin£ 
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(B4) 



(B5) 



(B6) 



2HD Model Realization: 



We discuss some details in context of the 2HD model example in Sec. VI. We set U(l)d charges as Qd[H{\ = 0, 
Qd[H 2 ] = Qd[Hd] = 1 for notational convenience. Then the gauge boson mass-squared is given by, with g z = 
g' / sm9w =g/ cos 9w, 

2 1 2 / 2 , 2\ 

m z« = - A 9z( v i+v 2 ), 



A At A x A\ - 

m z o = g d {v 2 + v d ) + 



e 1 2 1 
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*2 1 2 ! e 
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ff K + ^ 2 ), 



(B7) 



gzg'(v l +v 2 ). 



We assume m 2 z0 <§C m| which will be the case as long as (g^, egd, e 2 ) <§C g z and Vd is not exceedingly larger than 

d 

the electroweak scale. We define tan/? = v 2 /vi, tan /3d = v 2 /vd, and u 2 = v 2 + u| — (246 GeV) 2 . Then we have 

1 2 2 
- A 9zv , 

{A 2 ) 2 
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m 70 
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' Z" 



gi(vi + v 1 sin J j3 cos 2 0) = g\V 



' z" 



sin 2 /3 

~2 



sin f3 c 



(l-sin 2 /3sin 2 /3 d ), 



(B8) 



sin 2 f3 + e tan Qy/. 

9z 



10 



Gauge kinetic mixing e does not contribute to Zd mass but it affects the Z-Zd mixing angle £. 

(i) In the v-2 = limit (i.e. pure dark photon limit), the Zd mass is entirely from the Higgs singlet Hd and the 
Z-Zd mixing angle is provided entirely by e. We have 

m% ~ m| m% d ^ gjvj £ ~ e tan^^ (B9) 

which give 

M W m| 2 - £ta ? n ^™l V (BIO) 

\ —etanOwrriz m z +£ tan w^ / 



The mixing induced by the mass matrix cancels the effects because of field redefinition in Eq. ( A5) for the Zd induced 
neutral current coupling. 



(ii) In the e = limit {i.e. pure dark Z limit) 

j 
svn z B d 



ff> 2i ^(l-sin 2 /W/3 d ) 



■ n ■ n (BH) 

2g d . 2 o m z d sin /? sin p^ 

^ ~ sin p ~ — ■ . 

9z mz ^1 -sin 2 /3sin 2 /3 d 

Taking 1 — sin /3sin [3d — 1 is valid when |A 2 | <C m z am z o. In this limit 

to z ~ m z o m z d — m z° £ — e z (B12) 

with 

TYl y 

£z = -5 and 5 = sin 8 sin p^ . (B13) 

2. Higgs- Vector- Vector Couplings 



We assume no mixing among Higgs scalars and refer to the SM-like Higgs as H. From Eq. (Bl), we can get the 
relevant Higgs coupling to vector bosons. 

^scalar = -ChZzHZZ + CuZZ d HZZd + —Ch Z d Z d H ZdZd + ■ ■ ■ (B14) 

The Feynman rules for coupling of H to two vector bosons V\ and V2 are then given by ig^vCuVi v 2 ■ 
In the 2HD example, we get 

Chzz = C^ z cos/3(cos^ + £tan^ sm C) 2 

Crzz a = Chzz cos fi (cos£ + etan#;y sin£)(sin£ — e tan 9 w cost;) (B15) 

Cnz d z d = Chz Z cos /3(sin£ - e tan 6 W cos £) 2 



with C S H M ZZ = \g\v. 

The ratio of couplings is 



_ Cgzz d _ C H z d z d _ sing-gtanflyycosg 
Cf/zz Cff ZZd cos £ + e tan 9 W sin £ ' 



which, with small |£| <C 1 from Eq. (B8), yields 



2 9d „. 



showing that is not sensitive to e. 
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The relevant Higgs decay rates, for mz d *C to#, are given by 



T(H -> ZZ) 
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with couplings given in Eq. ( B15 ). Eq. ( B18 ) conveniently shows the effects of phase space and Z-Z d mixing in the 
Higgs decay rates. The ratio of Higgs decay rates in the Z d Z d and ZZ d channels is 



T(H -> z d z d ) e 2 



T{H -> ZZ d ) ' ' 2 m% d 
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(B19) 



where Eqs. ( Bll ) and ( B12 ) have been used in the second line. 



Appendix C: Zd Decay Width 



Using Eqs. ( |l0| ) and (11) in the text, we find that the partial decay width of Z d into the SM fermion pair // is 
given by, neglecting mf/mz d corrections Eqj, 



N c 



T(Z d -> //") ~ -^e 2 z g 2 z {$ f + g%) m Zd , 



48tt 



(CI) 



where g' v * = I3/ — 2Qf (sin 2 9\y — {z/sz) cos 9\y sin^vy) and gAf = — T 3 f. Here, Nc = 3 for quarks and Nc = 1 for 
leptons. 
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